Convective Conditions on MHD Casson Fluid in a Suspension of Dusty Particles with Darcy’s Porous
S. U. Mamatha1, C. S. K. Raju2*, D Madhavi Giridhar 3
1Department of Mathematics, Garden City University, Bengaluru, Karnataka 560049
2Department of Mathematics, GITAM University, Bangalore, Karnataka-562163
3Department of Mathematics, Garden City University, Bengaluru, Karnataka 577002
*Corresponding Author E-mail:mamathasupadhya@gmail.com1, sivaphd90@gmail.com2*, madhavi.g@gardencitycollege.edu
ABSTRACT:
The present study is to analyze Darcy-Forchheimer steady incompressible (MHD) non-Newtonian Casson fluid with suspension of dust particles over a stretching surface with exponentially decaying heat source. Convective boundary conditions are also considered. Similarity transformations were used to convert partial differential equations (PDEs) to a system of nonlinear ordinary differential equations (NODEs) which are solved numerically by employing Runge-Kutta with shooting technique. The effect of pertinent parameters on velocity and temperature profiles of both fluid and dust phase within the boundary layer has been studied by considering various values of controlling parameters. Additionally, skin friction coefficient and Nusselt number have been examined for various values of the governing parameters. It is found that non-uniform inertia coefficient of porous medium, porosity parameter and Casson fluid parameter decreases friction factor and Nusselt number.
KEYWORDS: Exponentially decaying heat source, Convective conditions, Dusty fluid, Casson fluid, Darcy-Forchheimer.
INTRODUCTION:
Over recent years, researchers are allured by the characteristics and performance of non-Newtonian fluids. Such motivation is due to complex physical behavior of non-Newtonian fluids. A non-Newtonian fluid finds numerous applications in pharmaceuticals, food products, crystal growth, fiber technology, chemical and material engineering. Hence several constitutive models for non-Newtonian fluids are available in the existing literature [1-4]. Non- Newtonian Casson fluid is a shear thinning fluid which exhibit yield shear stress. If suppose the yield stress is lesser than the applied shear stress than the liquid would start moving. For example tomato sauce, honey, human blood and fruit juices.
Casson fluid finds several applications in cancer homeo-therapy, fibrinogen, protein cells and red blood cells. Because of these applications researchers Raju et al.[5] studied heat and mass transfer of MHD Casson fluid on an exponential permeable stretching surface. Further Raju et al. [6] took up a comparative study of non-Newtonian Casson and Carreau fluid over a stretching surface and found that heat and mass transfer in Casson fluid is considerable high compared with Carreau fluid.
Phenomenon of suspended particles finds significant role in advanced processes of engineering and scientific field concerned with polymer technology, metallurgy, powder technology, power generation, atmospheric fallout, nuclear reactor cooling, combustion, waste water treatment, electrostatic precipitation of dust ,fluidization etc. Fundamental studies of suspended particles in the flow were initiated by Saffman [7]. Recently Krupa Lakshmi et al.[8] took up heat and mass transfer of a dusty liquid past a stretching sheet by considering thermal radiation in the energy equation. Gireesha et al. [8] studied influence of viscous dissipation on flow of dusty fluid over stretching sheet. Krishnamurthy et al. [9] took up theoretical study of suspended particle on slip flow and nanofluid in a porous Medium with the presence of nonlinear thermal radiation.
Extensive research has been undertaken by the researchers to study flow induced by a stretching surface in view of its large range of applications in chemical processing equipment, polymer extrusion, wire and fiber coating, glass fiber production, water pipes, irrigation channels etc. Due to its vast applications, initially, Sakiadis [11] presented the concept of two-dimensional boundary layer flows past a continuously stretching sheet with constant speed. Later on, several authors [12-22] extended the research over a stretching sheet for various fluid flow models over several geometries such as cone, Cylinder, Plate etc.
Heat and mass transportation from varied geometries embedded with porous media has several engineering and geophysical applications in geothermal reservoirs, thermal insulation, hydrology, agriculture, mineral processing, etc. The classical Darcy’s theory which has been used in most of the flow problems is suitable for lower velocity and smaller porosity. When inertial and boundary effects take place with higher flow rate Darcy’s theory is inadequate. Keeping this in view researchers [23-25] studied heat and mass transfer of fluid by considering Darcy- Forchheimer porous medium.
By studying the above literature, in this paper we proposed a mathematical model of boundary layer flow of Casson fluid with suspended uniform sized dust particles towards a stretching sheet in the presence of magnetic field, Darcy-Forchheimer porous medium, exponentially decaying space dependent internal heat source and convective boundary condition. The governing nonlinear partial differential equations are reduced into a set of non-linear ordinary differential equations by using local similarity transformation. These nonlinear ordinary differential equations then solved numerically by Runge-Kutta based shooting process for different values of parameters of interest. Influence of governing parameters on velocity, temperature, Skin friction coefficient and local Nusselt number are analyzed and discussed in detail with help of plotted graphs and tables.
Fig.1 Schematic representation of the flow diagram
2. MODELING:
A steady
two-dimensional, incompressible, magnetohydrodynamic boundary layer flow of Casson
fluid with suspended uniform size dust particles over a stretching sheet characterizing
Darcy-Forchheimerporous medium is developed. The stretching sheet is aligned
with the
axis and
the flow is confined to
.
The surface is stretched with the velocity
in
direction
with rate
.Magnetic
field of strength
is
imposed normal to the fluid flow as shown in Figure 1 By assuming
electrical conductivity of the fluid to be small we have neglected the induced
magnetic field. Convective surface temperature is characterized by
and heat
transfer coefficient by
.The
temperature
at the
surface of the sheet is considered to be more than the ambient fluid
temperature
. The
associated governing boundary layer equations are:
(1)
(2)
(3)
(4)
(5)
(6)
With boundary conditions
![]()
(7)
Here
- velocity
components along
direction
of the fluid and dust particle phase,
- velocity
components along
direction
of the fluid and dust particle phase,
are
density of the fluid, density of the dust particle phase, mass of the dust
particles, number density of the dust particle, kinematic viscosity, thermal
conductivity, electric conductivity, uniform magnetic field, specific heat of
the fluid and dust phase.
-permeability
of porous medium, A is power law index,
denotes
stokes drag coefficient.
-radius
of the dust particle.
is
the non-uniform inertia coefficient of porous medium, where
is the
drag coefficient [24] .
represents
temperature of the fluid, temperature of the dust particle phase, ambient fluid
temperature, convective fluid temperature, and convective heat transfer
coefficient.
and
represents
relaxation time of dust particles and thermal equilibrium time.
Coefficient
of the dimensionless space-dependent internal heat generation [26].
Following similarity transformation were used to convert PDEs into set of ODEs.
![]()
![]()
(8)
Continuity equations (1) and (3) are identically satisfied. Equations (2), (4)-(7) are transformed as follows.
(9)
(10)
(11)
(12)
The boundary conditions defined in equation (7) will be transformed to:
(13)
in the above expressions, represents magnetic
parameter, mass concentration of dust particles, fluid-particle interaction parameter,Casson
fluid parameter, porosity parameter, inertia coefficient, Prandtl number, heat
source parameter, Eckert number, specific heat ratio, fluid particle
interaction parameters,biot number.
(14)
Distinctive
measures of practical interest are skin friction coefficient
and local
Nusselt number
Which
are defined as ;

Where
is the local
Reynold’snumber.
Fig.2 The effect
of
on
velocity field Fig.3 The
effect of
on
temperature field
Fig.4 The effect
o
on
temperature field Fig.5 The
effect o
on
velocity field
Fig.6 The effect
o
on
velocity field Fig.7 The effect o
on
velocity field
Fig.8 The effect
o
on
temperature field Fig.9 The effect of
on
temperature field
Fig.10 The effect
o
on
temperature field
Table 1
Comparison of the results for local Nusselt number with
.
|
Pr |
Ishak [17 ] |
El-Aziz [8 ] |
K.L.Krupa Lakshmi et.al[2] |
Present Study |
|
0.72 |
0.8086 |
0.80873 |
0.808630 |
0.8075 |
|
1.0 |
1.0000 |
1.0000 |
1.0000 |
1.0000 |
|
3.0 |
1.9237 |
1.92368 |
1.92367 |
1.9146 |
|
10.0 |
3.7207 |
3.7207 |
3.72067 |
3.6516 |
|
100 |
12.2941 |
12.2941 |
12.294087 |
12.2936 |
Table 2 Numerical
values of Skin friction
and
Nusselt number
for
different values of
.
|
|
|
|
|
|
|
|
|
|
0.5 |
|
|
|
|
|
-0.442845 |
-0.344656 |
|
1 |
|
|
|
|
|
-0.521758 |
-0.350769 |
|
1.5 |
|
|
|
|
|
-0.591849 |
-0.356023 |
|
|
0.1 |
|
|
|
|
-0.442845 |
-0.209401 |
|
|
0.3 |
|
|
|
|
-0.442845 |
-0.440393 |
|
|
0.5 |
|
|
|
|
-0.442848 |
-0.566926 |
|
|
|
0.1 |
|
|
|
-0.430626 |
-0.343955 |
|
|
|
0.3 |
|
|
|
-0.454752 |
-0.345334 |
|
|
|
0.5 |
|
|
|
-0.477702 |
-0.346621 |
|
|
|
|
0.1 |
|
|
-0.425760 |
-0.343328 |
|
|
|
|
0.3 |
|
|
-0.459462 |
-0.345950 |
|
|
|
|
0.5 |
|
|
-0.491397 |
-0.348431 |
|
|
|
|
|
0.1 |
|
-0.331109 |
-0.338235 |
|
|
|
|
|
0.3 |
|
-0.517415 |
-0.348701 |
|
|
|
|
|
0.5 |
|
-0.616068 |
-0.353715 |
|
|
|
|
|
|
0.1 |
-0.442845 |
-0.355359 |
|
|
|
|
|
|
0.3 |
-0.442845 |
-0.335014 |
|
|
|
|
|
|
0.5 |
-0.442845 |
-0.318266 |
RESULTS AND DISCUSSION:
For analyzing
the approximate solutions of velocity
as well as
temperature
fields,
the non-linear ordinary differential Eqs. (9)-(12) with reference to the
boundary conditions (13) are solved numerically with the assistance of
Runge-Kutta and Newton’s methods. For numerical solutions we have considered
the values of non- dimensional parameters as![]()
![]()
these
values are taken constant in this study besides the varied parameters as
mentioned in the respective figure. In Figure 2and 3 the influence of
magnetic field parameter
on
velocity
and
temperature profiles
of
fluid and dust phases can be observed. As
accelerates,
The Lorentz force which opposes the flow regime also improves hence there is
decrement in the
and
Lorentz force generates friction on the flow, which results more heat energy
causing increase in the temperature distribution of the flow.
Figure 4
demonstrates the influence of convective heating parameter
on the
temperature distribution of fluid and dust phase. It is noticed that increasing
values of
improves
thermal boundary layer
.
Physically convective heat exchange at the surface of the sheet helps to
improve the thermal boundary layer thickness. Figure 5 presents larger
values of non –uniform inertia coefficient
declines
the velocity profiles of both fluid and dust phase. Figure 6 and 7
depicts the variation in velocity and temperature profiles for larger values of
porosity parameter
.
It is observed that increment in
reduces
velocity field and improves temperature field for both fluid and dust phase. In
Figure 8 and 9 it is noticed that increment in Casson fluid parameter
decreases
velocity boundary layer and improves thermal boundary layer for both fluid and
dust phase. Figure 10 demonstrates the influence of exponentially
decaying internal heat generation parameter
over
temperature profiles of fluid and dust phase. It is observed that increase in
enhances
heat generation which causes the temperature profiles of fluid and dust phase
to enhance.
Table 1 Shows the comparative study of thepresent result with the published results. Present results agree with the published results.
Table 2
display the effect of non-dimensional parameters
on skin
friction coefficient
and
local Nusselt number
.
It is noticed that rise in magnetic field parameter minimize the friction
coefficient along with Nusselt number. Increment in Biot number decrease
Nusselt number. But no variation is observed in friction factor. Rising values
of inertia coefficient, porosity parameter and Casson fluid parameter decreases
skin friction and Nusselt number. Increment in exponentially decaying internal
heat generation parameter enhances heat transfer rate and friction factor
remains same.
4. CONCLUSION:
Darcy –Forchheimer two- dimensional flow of non-Newtonian Casson fluid with suspended uniform size dust particles over a stretching sheet with exponentially decaying heat generation and convective boundary condition has been discussed. The major findings of present study are as follows.
· The magnetic interaction parameter motivates friction factor coefficients and heat transfer rate.
· Increment in Biot number improves the thermal boundary layer of both phases. But retards heat transfer rate.
· Heat source parameter declines thermal boundary layer but improves heat transfer rate.
· Non-uniform inertia coefficient, Porosity parameter, Casson fluid parameter reduces friction factor and heat transfer rate but improves thermal boundary layer in both the fluid and dust phase
5. REFERENCES:
1. Raju CSK, Sandeep N, Jayachandra Babu M, Sugunamma V. Dual solutions for three-dimensional MHD flow of a nanofluid over a nonlinearly permeable stretching sheet. Alexandria Eng. J. 2016; 55:151–162. doi:10.1016/j.aej.2015.12.017.
2. Mamatha SU, Mahesha, Raju CSK.Multiple Slips on drodynamic Carreau Dustynano Fluid Over a Stretched Surface with Cattaneo-Christov Heat Flux. J. Nanofluids. 2017;6: 1074–1081
3. Mamatha S U, Raju CSK, Madhavi G , Mahesha. Unsteady 3D MHD Carreau and Casson Fluids over a Stretching Sheet with Non-Uniform Heat Source/Sink.IISTE. 2017;52: 10-23
4. Raju CSK, Sandeep N, Sulochana N, Babu M.J. Dual Solutions of MHD Boundary Layer Flow past an Exponentially Stretching Sheet with Non-Uniform Heat Source / Sink.J. Appl. Fluid Mech. 2016; 9(2). 555-563 .doi:10.18869/acadpub.jafm.68.225.24784.
5. Raju CSK, Sandeep N, Sugunamma V, Babu M.J, Reddy JVR, Heat and mass transfer in magnetohydrodynamic Casson fluid over an exponentially permeable stretching surface. Eng. Sci. Technol. an Int. J. 2016 ;19 :45–52. doi:10.1016/j.jestch.2015.05.010.
6. Raju CSK, Sandeep N, Unsteady three-dimensional flow of Casson –Carreau fluids past a stretching surface.AlexandriaEngineeringJournal.2016;55(2):1115-1126 http://dx.doi.org/10.1016/j.aej.2016.03.023.
7. Saffman PG. On the stability of a laminar flow of a dusty gas. J. Fluid Mechanics.1962;131:20–128.
8. Krupa Lakshmi KL, Gireesha BJ, Rama S R Gorla, Mahantesh B.Two-phase Boundary Layer Flow , Heat and Mass Transfer of a Dusty Liquid past a Stretching Sheet with Thermal Radiation. Int J. Industrial Mathematics .2016;8(3): 14 pages.
9. Gireesha BJ, Roopa PG,.Bagewadi CS.Effect of viscous dissipation and heat source on flow and heat transfer of dusty fluid over unsteady stretching sheet. Appl Math Mech. 2012; 33: 1001-1014.
10. Krishnamurthy MR., Gireesha BJ, Gorla RSR, Prasannakumara B.C. Suspended Particle Effect on Slip Flow and Melting Heat Transfer of Nanofluid Over a Stretching Sheet Embedded in a Porous Medium in the Presence of Nonlinear Thermal Radiation. J. Nanofluids. 2016; 5:502–510.
11. Sakiadis BC, Boundary Layer Behavior on Continuous Solid Surfaces:II Boundary Layyer on a Continuous Flat Surface. AIChE Journal. 1961; 7(2):221-225.
12. Crane LJ. Flow past a stretching plate. ZAMP.1970; 21645-647.
13. Raju CSK,.Sandeep N, Sulochana C, Sugunamma V, JayachandraBabu M , Radiation, inclined magnetic field and cross-diffusion effects on flow over a stretching surface. Journal of the Nigerian Mathematical Society. 2015); 34(2):(169-180.
14. Raju CSK, Sandeep N,. Jayachandra Babu M, Sugunamma V, Dual solutions for three-dimensional MHD flow of a nanofluid over a nonlinearly permeable stretching sheet. Alexandria Eng. J. 2016;55: 151–162. doi:10.1016/j.aej.2015.12.017.
15. Jayachandra BabuaM, Raju CSK,Sandeep N, SulochanaV, Sugunamma C,. Radiation, inclined magnetic field and cross-diffusion effects on flow over a stretching surface. Journal of the Nigerian Mathematical Society. 2015;34(2):169-80.
16. SandeepN, Sulochana C, Raju CSK, Jayachandra Babu M, Sugunamma V.Unsteady boundary layer flow of thermophoretic MHD nanofluid past a stretching sheet with space and time dependent internal heat source/sink.Appl. Math.2015;10(1):312-327
17. Sohail Nadeem a, Rizwan Ul Haq a, Zafar Hayat Khan,Heat transfer analysis of water-based nanofluid overan exponentially stretching sheet.Alexandria Engineering Journal.2014; 53: 219–224
18. Reddy JV, Sugunamma V, Sandeep N, Raju CSK.Induced Magneticfield Effect on Stagnation-Point Flow of Magneto-Nanofluids Towards a Stretching Sheet.Advanced Science, Engineering and Medicine.2015; 7(11).968-974
19. Ishak A, Nazar R, Pop I. Heat transfer over a stretching surface with variable heat flux in micro polar fluids. Phys. Lett. A. 2008; 372(5): 559–61.
20. Aziz MA, Mixed convection flow of a micropolar fluid from an unsteady stretching surface with viscous dissipation. Journal of the Egyptian Mathematical Society. 2013 ;21:385-394
21. Raju CSK, Sandeep N, Babu MJ, Reddy JVR. Stagnation-Point Flow of a Ferrofluid Towards a Stretching Sheet.Journal of Nanofluids. 2016;5: 1–8. doi:10.1166/jon.2016.1209.
22. Mamatha SU, Mahesha, Raju CSK .Cattaneo-Christov on heat and mass transfer of unsteady Eyring Powell dusty nanofluid over sheet with heat and mass flux conditions. Informatics Med. Unlocked. 2017;9:76-85doi:10.1016/j.imu.2017.06.001.
23. Jumah RY, Fawzi A, Abu-al-rub F, Darcy- Forchheimer mixed convection heat and mass transfer in fluid saturated porous media. International Journal of Numerical Methods for Heat and Fluid Flow. 2001;11(6): 600 – 618
24. Hayat T, Haider F, Muhammad T, Alsaedi A. Three-dimensional rotating flow of carbon nanotubes with Darcy-Forchheimer porous medium. Plos/one. 2017 ;1–18.
25. Seddeek MA, Influence of viscous dissipation and thermophoresis on Darcy – Forchheimer mixed convection in a fluid saturated porous media.Journal of Colloid and Interface Science. 2006;293 :137–142. doi:10.1016/j.jcis.2005.06.039.
26. Animasaun Isaac Lare. Casson fluid flow with variable viscosity and thermal conductivity along exponentially stretching sheet embedded in a thermally stratified medium with exponentially heat generation. Journal of Heat and Mass Transfer Research.2015;2:63-78
Received on 11.06.2017 Modified on 12.07.2017
Accepted on 18.09.2017 ©A&V Publications All right reserved
Research J. Science and Tech. 2017; 9(3):453-460.
DOI: 10.5958/2349-2988.2017.00079.1